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SUMMARY 

Under the assumption that a wing, body, or wing-body combi- 
nation is slender or flying at near sonic velocity j expressions are 
given which permit the calculation of pressure in the immediate 
vicinity of the configuration. The distarhance field, in both 
svisonic and super sonic dight, is shown to consist of two-dimen- 
sional disturbance fields extending laterally and a loTigihidijial 
field that depends on the streamwise ^owth of cross-sectiorial 
area. A discussion is also given of couplings, between lifting 
and thickness eflecis, ihai necessarily arise as a result of the 
quadraiic dependence of pressure on the induced velocity 
components. 

INTRODUCTION 

This paper is concerned with the prediction of pressure 
distribution on or in the immediate vicinity of a wing, body, 
or wing-body combination under conditions in which the 
geometric configuration is slender in the flight direction or 
is flying at near sonic velocity. The material to be presented 
is thus associated with the rather extensive group of re- 
sults that belong to what is usually referred to as slender- 
wing theory. The basic assumptions and methods can be 
found in publications by Munk, R. T. Jones, and Ward 
(refs. 1, 2, and 3) ^ and a discussion of the applicability of 
the methods to the prediction of loading on slender wings at 
sonic flight speeds has been given in reference 5. In refer- 
ence 2, attention was directed toward the calculation of load 
distributions over wings in subsonic and supersonic flight 
and reference 3 was devoted to the consideration of super- 
sonic flight velocities. It is therefore of interest to investi- 
gate further the effects attributable to thickness on wings 
and wdng-body combinations at both subsonic and super- 
sonic flight speeds. Such investigations lead to valid ap- 
proximations of interference effects and also indicate the way 
in which thickness and lifting effects can produce couplings 
in the calculations of pressures induced in the flow fleld. 

ANALYSIS 

It is proposed to take the basic solutions of the linearized 
partial differential equations governing three-dimensional 
compressible flow and to obtain a simplification of the ex- 
pressions by restricting attention to the induced field in the 
immediate vicinity of slender airplanes or missiles. These 

i Supersedes NAOA TN 2900, *‘The Oalcnlatlon of Pressure on Slender Airplanes In Sub- 
sonJo and Supersonio Flow,^' by Max, A. Heaslet and Harvard Lomax, 1963,* 

* Reference should also be made to the extensions of slender-wlng theory by Adams and 
Sears (ret 4). 


simplified expressions contain solutions used previously to 
study the forces and moments on lifting wings and bodies. 
In addition, however, they can be used to evaluate the 
first-order thickness effects on the pressure in the vicinity of 
the wing and body. 

Consider, first, the construction of a weakly disturbed 
flow field. Let a uniform stream flow in the direction of 
the positive x axis of a Cartesian coordinate system, as in 
figure 1. Immerse in the stream, which has a velocity Z7o 




Figure 1. — Orientation of coordinate system and perturbation 

velocities. 

and a Mach number Mo, a slender wing-body shape the 
surface of which is inclined at a small angle to the free- 
stream direction. This angle of inclination must be small 
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enougli so that nearly everywhere in the fluid the magnitude 
of the perturbation velocity vector divided by the speed of 
the free stream is much less than one; that is, 

(la) 

Moreover, large supersonic Mach numbers are to be avoided 
and as a measure of this condition the inequality 

(lb) 

is imposed. 

Consider, next, the linearized partial differential equation 
governing weakly disturbed isentropic fluid flow. In terms 
of the perturbation velocity potential (p{Xy y, 2 ), the lowest 
order approximation consistent with inequalities (la) and 
(lb) is ' ’ 

(1 , ( 2 ) 

where the subscripts denote partial differentiation with 
respect to the indicated variable. 

Consider, finally, the expression for the pressure coefficient 
that is again consistent to the lowest order with inequalities 
(la) and (lb). By expanding the pressure-velocity relation 
for steady isentropic flow and n^lecting higher-order terms, 
one finds 


p _ p— Po 2u (1— 

1 u, ^ w 

where p and p are pressure and density, respectively, and 
the subscript 0 refers to conditions in the free stream. It 
follows from inequalities (la) and (lb) that pressure co- 
efficient can be expressed in the form 




2u 

Uo W 


(3) 


Equation (3) is the simplest general expression for pressure 
coefficient that is stiU entirely consistent with the assump- 
tions basic to the development of equation (2). 

Special solutions applying to problems of the class indi- 
cated can be obtained by appropriate simplification of general 
solutions to equation (2). Such a procedure will be discussed 
in the next section. The pressure coefficient is then deter- 
mined by substituting these results into equation (3). The 
simplifications that can be made in evaluating the pressure 
on the surface of the airplane will also be discussed. 

THE REDUCED SOLUTIONS 

Subsonic. — As it apphes to subsonic flow, equation (2) 
can be written in its normalized form as 


+ <Pxx=0 


(4) 


A well-known solution to equation (4), resulting from an 
application of Greenes 'theorem, is given by the expression 


tJsIQn' 


where dSi is the element of surface area on the airplane or its 
vortex w^e, r equals V(^ ^ 1 )^} aiid 5/5^' is the 
derivative normal to the surface Si. When this solution is 
apphed to boundary-value problems for slender configurations 
it can be simplified considerably. 

Eor example, when the airplane shape is slender it is 
justifiable to iutroduce simplifications in the form of the 
derivative b/dn' and the differential aread^'i. The operator 
b/bn^ can be expressed as 


b , b , b 

where %, and % are the direction cosiues between a normal 
to the surface ^L^iid the x, y, and z axes, respectively. The 
differential area dSi can be expressed as 

dsidxi 


where dsi is a differential length along the surface in a 2/2 
plane. If the shape is slender, Ui is small and can be neglected 
relative to either unity or 

By means of these simplifications, equation (6) can be 
approximated by the expression (from now on, the configura- 
tion will be considered to he along. the positive Xi axis Avith its 
foremost part iu the a:i=0 plane) 




dsi 




( 6 ) 


where bfbn represents n^fby-\-n^lbZj ’the normal derivative 
to a section in the yz plane, and s is the cmve boimdiug this 
section. 

If the wing-body configuration is slender, the ratio 
\r[{x—x^Y is small over almost aU of its surface and vortex 
wake provided the point x^y^zis on or in the vicinity of these 
surfaces. This imphes the approximation 

■yj{x—XiY+T^^\x—Xi\ (7) 

can be used -to simplify further equation (6). However, 
since, in the limiting case of r=0, equation (6) is a divergent 
integral, it is necessary to introduce this approximation 
with some care. 

First let us consider in equation (6) only the portion of the 
integral multiplying (pbfbn. Designating this by (piix^y^z) 
one can readily show 


<pi{x,y,z) 


AA r ip(xiy8^{x—xi)(^rlbn) 

47t bxjo V, 


dsi (8a) 


The analysis of equation (4) can be interpreted as applying 
to the condition Afo=0 but one can extend the solutions 
throughout the subsonic Mach number range by applying 
the Prandtl-Glau^t rule. 


which, with the approximation given by equation (7), re- 
duce to .. 

fM.y.c) ” l£lf ^ 
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Since this can be written 


^ (J‘ fc-X *■)// S *■ 

it simplifies to 

<Pi(x,y,z ) « ^ 


(9) 


The other term in equation (6) requires more attention. 
Designating the velocity potential induced by this term as 
<ptj we can write 


<pi(p:,y,x)=- 


J_ A 

47T bx 



b<p X—Xi 

bn 


r 


( 10 ) 


The logarithm in the integrand separates into two parts and 
if the term containing ln{l!r) is further simplified, in the 
maimer used to derive equation (9) from equation (8) , one has 




_L A f” 

4tr dxjo \x—Xi\ 
■yj{x—x{f+r^\ dsT. 


dxi^^^ln{\x-x,\+ 


( 11 ) 


To the second term in equation (11) we now apply the 
mean-value theorem. First divide s, the curve bounding the 
airplane^s normal cross-sectional area, mto so that 

bip/bn has the same sign everywhere along each Sj. This Has 
a clear physical interpretation, for, since <p is the perturbation 
velocity potential, to the order of our approximations 




( 12 ) 


where Sj{xi) is the part of the cross-sectional area with its 
only exterior boundary along the arc S/. (Since we are con- 
cerned only with the rate of change of Sj(x)j the internal 
boundary is immaterial.) Hence, if bip/bn is everywhere 
positive along 5 /, the surface is everywhere expanding there; 
and, conversely, if b<pjbn is negative, the surface is con- 
tracting. 

Using these definitions, we can apply the mean-value 
theorem to the second term in equation (11) and write 


J^|2z»[lx-Xxl+V^=^l?+^ dsi=Uo'^ln[\x-xi\+ 

■yj{x—x{)^+rj‘\ S'j(Xi) 

(13) 

where r,=-yj{y—yjy+(z—zjy and yj, % is a point on the 
surface of the configuration in the a?i plane somewhere along 
Combining the above resul ts, and a pplying the Prandtl- 
Glauert transformations 03=V|Afo*“-ll, V^PVi 
one can now approximate equation (6), when it is applied to 
the flow field in the vicinity oj slender wing-body crnifiguraiumSy 
and further, when it is applied to conJigvrcUions for which 
S'{x) is continuous, and for which S"j(x) exists, by the equation 




(14) 


Notice that the limitation of equation (14) to the vicinity 
of slender shapes has, as yet, entered only in the approxima- 
tion of equation (8a) by equation (8b), that is, in approxi- 
mating the effect of the (pbfbn{lnr) (or doublet) term as given 
in equation (14). If equation (14) is applied to the study of 
thin, nonlifting, uncambered wings, therefore, it is only 
limited by the assumption that (x) is continuous. Further, 
if such wings are slender, the position of the point yj, Zj does 
not deviate far from zero. Hence, we can chose for the 
value ro, where Then, since ^.S7 (x)=S" (x) . 

integrating equation (14) by parts yields 

%S'ix) Inro 

(15) 

which is the result presented by Keune in reference 6. 

If we continue to study the flow in the vicinity of general 
slender shapes, however, equation (14) can be further sim- 
plified by applying the approximation given in expression 
(7) . In the first place, if S"(x) is continuous (and, therefore, 
vanishes at x=0 and Z, Z being the total airplane length), the 
equation for the potential can be approximated everywhere 
in the vicinity of the configuration by the simple expression 

p(x, y, z)^<p3(x; y, 2 )— S”(xi) In 

(16) 

where S{pi) is the total cross-sectional area in a plane normal 
to the free stream and where ® 


/ , 1 ^ UoT S\xi 


(xi)dxi 




< 

On the other hand, if S"{x) has a discontinuity, the approx- 
imation given by equation (16) yields a logarithmic discon- 
tinuity in the value of bipfbx, that is (see eq. (3)), in the 
pressure coeflB.cient. Tins discontiauity is spurious (it does 
not exist in solutions given by exact linearized theory) and 
it can be avoided by modifying slightly the simplifying 
procedure followed above. 

Suppose, for example, S"{x) has a finite discontinuity 
of magnitude AjS}^(Zi) along the arc Si in the plane x=h, 
so that 


S"{x)=\ (18) 

is:'{x)+ASy(liy, h<x 


* Notice that p, *) Is a solatlon to Laplace's equation in two (the v and z) dimensions. 
The X dimension d^ not appear explicitly In this part of the oomplete solution for p, i) 
but enters as a parameter when w is adapted to particular boundary conditions. 

i Any integrable singularity in S^'(r) could be treated but the aualsrsls was restricted for the 
sake of simplicity. 
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where S”(x) is everywhere continuous and AS'/(h) is a 
simple step function. Then equation (14) reduces to 


<p(x, y, z) « y, 2 )— In dxi— 

47t Xi| L P 


^dxi 

(19) 


Further, when x lies between 0 and Z, equation (19) can be 
written (since Sc'(xi) — —ASy(li) for xC>F) 


<p(x,y,z)^<Pi(x;y,z)—^j^ Sl(xi) In ^ dxi— , 


UoAS’,(li) 

4v 


{(Z-a;) 


( 20 ) 


The physical significance of equations (16) and (20) is 
clear. If the second derivative of the area distribution is 
continuous, the three-dimensional velocity field induced by 
airplanes that are slender is approximated in the vicinity of 
their surfaces, by 

(1) a velocity field that is independent of the Mach num- 
ber and satisfies the two-dimensional Laplace equation 
and boimdary conditions in transverse planes 

(2) a longitudinal field that depends on the Mach number 
and streamwise variation of cross-sectional area and 
is independent of y and z. 

If the second derivative of the cross-sectional area is dis- 
continuous, part (2) of the above interpretation should be 
modified to apply only to the continuous portion of S"{x) 
and additional terms, depending on the number and location 
of discontinuities in jS"(a;), are included. 

Supersonic. — ^In the case of supersonic flow, the normalized 
form of equation (2) becomes 


^xx— <pzx=0 ( 21 ) 

An analysis based on equation (21) applies specifically to the 
condition Afo=V2 but these results can be extended through- 
out the supersonic Mach number range by applyiog the 
Prandtl-Glauert rule. Volterra’s solution to equation (21), 
(see, e. g., ref. 7, p. 190) which is analogous to the subsonic 
form given in equation (5), is expressible as 

■pM — 

( 22 ) 

where, as in the subsonic case, dSi is an element of surface 
area on the airplane or its vortex sheet and r equals 
V(y — 2/i)^+(2— 2i)^* ^ distinction to the subsonic solu- 
tions, the area of int^ation is now the portion of the 
airplane and its vortex wake within the forecone from the 
point Xjy,z and b/dv is the derivative along the conormal * 
rather than the normal. 


^ The conormal, for tbo normalized form of the differential equation, is the vector that 
results from chang in g the algn of the x component of the n ormal . 


If the conormal and dififerential area are expressed in 
terms of the direction cosines and the application of equation 
(22) is limited to slender configmations, rti can again bo 
neglected relative to unity or Furthermore, the 

approximation, similar to expression (7) for the subsonic 
case, 

V(a^— Xi)*— a^il 

is implied. Under these conditions, the potential in the 
vicinity of slender shapes flying at supersonic speeds can be 
approximated by the equation 

<p(x,y, z ) « 4 :) 

gs SY(x^ Z7.[£z£i±V^lBVj^^ 

(23) 


where, again, S'(x) must be continuous. Equation (23) 
diBfers from equation (14) only by a factor of 2 in the second 
term and the extent of the Xi integration. In the supersonic 
case the Xi integration is carried only to x—firjj or to x when 
Tj can be neglected, smce the original integration area r in- 
cluded only the points in the forecone from Xjy^z (these two 
differences were compensating in the derivation of the first 
integral term). The second term in equation (23) further 
simplifies in a manner analogous to that used for the simpli- 
fication of equation (14). For example, the expression for 
the perturbation potential near a slender configuration hav- 
ing a discontinuity in S^'(x) along the arc Sj in the plane 
x=li is given by 


v(x,y, z ) « pa(x;y, z) — ^ S',' (xi) In ^ (Zxi+ 

Uo^S’/(h)f . x—h—prj\( j 

(Zi-x) In | 

(24) 


where S'J (x) and AjS7(Zi) are defined by equation (18) and 
^(x,*y, 2 ) is the two-dimensional solution to Laplace^s equa- 
tion defined by equation (17). 

The physical significance of equation (24) is analogous to 
that for equation (20). Before proceeding to the next section, 
however, two observations regarding these solutions are 
worth mentioning. In the first place, notice that if discon- 
tinuities in S'/ {x) occur on the x axis, as could be the case, 
for example, at the nose of a pointed body of revolution or at 
the apex of a triangular wing, and if one is interested in 
evaluating ^(a;,y, 2 r) only on the objectis surface, the value of 
Ti for such discontinuities would be zero and, if there were no 
other discontinuities, equations (20) and (24) would be 
correct if the AjS 7 '(Zi) were set equal to zero and S/{^\) was 
written simply S^^{x/). In the second place, for the super- 
sonic case, only the discontinuities between the nose and the 
plane x^ at which the induced velocities are being calculated, 
affect the flow there. Hence, unless one is interested in the 
flow field behind the configuration, any discontinuity at 
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x^l can always be excluded from equation (24). Taking 
these observations into account simplifies the treatment of 
many body-of-revolution and supersonic-wing problems. 

THE BEFERENCE COORDINATE SYSTEMS 

Equation ( 2 ) was developed specifically for the case in 
which the tmdisturbed stream at infinity is parallel to the x 
axis. A coordinate system so oriented is usually referred to 
os the wind-axes system. (See fig. 2 .) When the configura- 



(a) Body axes. 

(b) Wind axes. 

Figure 2. — ^Illustration of wind- and body-axes sj^stems. 


tion is tilted with respect to the free-stream vector, however, 
it is often easier to study the boundary-value problem with 
axis placed along the center line of the fuselage. Such a 
coordinate system is usually referred to as the body axes. 

Obviously the wind and body axes differ significantly 
only by rotations about the y and z axes. When Mq is 
zero, equation ( 2 ) is invariant to such a rotation, but for 
values of Mo greater than zero this is no longer true. How- 
ever, when Mo is greater than zero, equation ( 2 ) represents 
the governing differential equation only to a certain order, 
and, if the magnitude of the rotation is similar to that of the 
parameters by which the equation is ordered, it is, in this 
sense, still invariant to rotations about aU three axes for 
both subsonic and supersonic Mach numbers. Thus equa- 
tion ( 2 ) is to the lowest order the governing partial differ- 
ential equation for both wind and body axes, provided the 
airplane is slender and the angles of attack and sideslip 
are small. 

Although the partial differential equation is invariant 
with respect to a small rotation of the coordinate system, 
the boundary conditions and expression for the pressure 
coeflS.cient in terms of the perturbation velocities are not. 


The following will contain a discussion of the boimdary con- 
ditions and the pressure and loading coefficients with 
reference to a body axes system. 

THE BOUNDARY CONDITIONS 

The boundary conditions require that the gradient of the 
total velocity potential evaluated infinitely far from the 
aircraft be consistent with a uniform free stream there (the 
direction of which depends on the orientation of the coordi- 
nate S 3 ^tem) and when evaluated normal to and on the 
surface of the airplane itaeff be zero. Let 4 >(z,y, 0 ) denote 
total velocity potential, ??(a;,y, 2 !) perturbation velocity poten- 
tial, and refer the analysis to body axes in a free stream. 
If the orientation of the free-stream velocity vector to the 
system of axes is fixed by the angles a and 7 as shown in 
figure 2 , one can write 

$(x, 2 /, 0 ) = TJq{x cos a cos y+y sin y+z cos y sin 
such that on the aircraft surface 
C/o(% cos a cos sin y^n^ sin a cos y) +ni^+7i2^3y+ 

Til, ^ 2 , and nz again being the direction cosines of a normal 
to the airplane surface with respect to the x, y, and z axes, 
respectively. By the assumptions basic to the present 
theory, the latter equation reduces to 

Uo(ni+7i2y+7iza)+‘~ <p{x,y,z)=0 (26) 

where, as before, n is the normal to the curve bounding a 
cross section in the yz plane. 

Equation (26), which applies to arbitrary slender shapes, 
can be simplified for many specific problems. Consider 
now three types of configurations that lead to such simpli- 
fications: first, a surface, such as a wing, which deviates 
'‘only slightly from a plane; second, a surface which forms a 
body of revolution; and, third, a surface which is a com- 
bination of the above two. 

Planar systems. — ^Let h{x,y) be the "distance a surface 
deviates from the z =0 plane, and s be local semi-span. 
(See fig. 3 (a).) Assume that aj{d8ldx)<i<Cl holds; then 
furthermore, if the inequality (dA/5x)/(cfo/(&)<^<^l is satis- 
fied, it is consistent with the previous approximations to 
neglect the y component of the normal along the wing sur- 
face and to project the velocity vector represented by the 
resulting vertical derivative to the upper or lower surface 
of the 2=0 plane. In this way equation (25) becomes 

Z7o(7li+7^3a)+Tl3 ^—0 

and, since 711 / 713 = —bA/ba;, the boundary conditions for planar 
problems ® are expressed by the equation 

df)...— 


« Certain planar systems, such as the cruciform \ving, require more than one plane but the 
concepts are essentlaUy the same as those presented here. 
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(b) Bodies of revolution. 

(c) Wing-body combinations. 

Figuhb 3. — Types of configurations for which boundary-value 
problems are discussed. 

Bodies of revolution. — ^Let R be the radius of a body of 
revolution. (See fig. 3(b).) If 0 is measured from the z axis 
in the yz plane and 7 is set equal to zero, the relations 


• The rather' obvious extension of the above concepts is to 
apply equation (26) over the winged portion of the con- 
figmution and equation (27) over the body. It is then 
necessary, however, to consider the relative magnitudes of the 
terms bhfdx, ds/dx, dtjdx and dRIdx, since they appear in the 
I solutions in various combinations: If the winged portion 

! is to be treated as a planar problem, the magnitude of 
bhfdx must be small enough to be neglected in comparison 
to the leading- and trailing-edge slopes, ds/dx and dt/dx. 
But this does not imply that bh/bx can be neglected in com- 
parison to dR/dx or that dR/dx can be neglected in com- 
parison to either ds/dx or dt/dx. The latter approximations 
will not, in general, be made. 

THE PRESSURE COEFFICIENT 

The expression for the pressure coefficient given by 
equation (3) is written in terms of velocity components 
that are referred to the wind axes. Its re-expression in 
terms of velocities referred to the body axes is readily 
determined. For the orientation shown in figure 2 (a), the 
equation becomes 

Op^ «v?») (28> 

Equation (28) can be used, in general, to evaluate the pres- 
sure in a perturbation velocity field that is refennd to the 
body axes. If the interest is limited to the pressure on the 
surface of the aircraft, however, certain simplifications can 
be made. For example, consider the configm’ation illus- 
trated in figure 3 (c) consisting of a sweptback wing moxmtod 
on a body of revolution. For simplicity, let 7=0. Applying 
the boundary conditions given by equations (26) and (27), 
one can show that on the surface of the wing 

and on the surface of the body 

«)!+“■-(§)’ » 

These solutions can be simplified further by considering the 
detailed nature of the perturbation velocity field induced by 
shapes such as that shown in figure 3 (c). For example, if 
S^'(x) is continuous, the results given by equations (20) and 
(24) can be expressed in the form 


^dR/dx dR cos 6 

7ll= - . » — j-y nz= —r — : 

■yfi+WJd^ ^ 


« cozd 


together with equation (25), give for the normal derivative 
on the surface of the body 




C7o cos e 

B ^ 


(27) 


which is the simplified expression of the boundaiy conditions 
pertaining to bodies of revolution. 

Interference effects. — Consider, finally, surfaces that are 
a combination of the above two as, for example, the one 
shown in figure 3 (c). 


<p{^yyy^)=<P 2 {x;y,z) +A(x) (30> 

where the expression for A(x) depends on whether the speed 
is subsonic or supersonic. Further, for the particular con- 
figurations being considered, the expression for <p%{x^j^z) can 
be written in the general form 

( 31 )' 

since the dependencgr on x can enter only through the bound- 
ary conditions •which, in turn, axe specified by the body 
radius B(,x). the -wing thickness, h{x,y), and the lateral 







THE CALCTJIATION OF PRESSURE ON SLENDER AIRPLANES IN SUBSONIC AND SUPERSONIC FLOW 


651 


distances from the center line to the trailing edge and leading 
edge, t{x) and s{x)^ respectively. The term a<pa will he re- 
ferred to as the potential due to angle of attack, since it 
vanishes when the angle of attack vanishes and increases 
linearly with increasing a; the term 

will be referred to as the potential due to thickness, since it 
exists when the angle of attack is zero, does not change with 
angle-of-attack change, and vanishes when the thicknesses 
of the wing and body do not vary with x,' 

By breaking ^ down into its component parts as in 
equation (31), it has been ordered in that the magnitudes of 
the terms on the right-hand side of equation (31) are con- 
trolled by the coefficients of the ^^s, and the derivatives of 
<pa, <Pbj and <pe with respect to 5, t, 2?, and z can aU be con- 
sidered equal. Since a and bh/dx are negligible relative to 
(It/dx and (h/dx (as was pointed out in the discussion of the 
boundary conditions for interference problems), equations 
(29a) and (29b) can be written: ' 

on the surface of the wing 


^ ( 2 b(p I 1 diR> ^ b(pa b(pc 

\Uo dx L by by 


+ 


by by~^dx\byjj],. 
and on the surface of the body 


(32a) 


a 


1 dE 


_ (2 b(p /dBV 

\Uo bx^\dxj dx 
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+ 




(32b) 


O, 


2 

' L UobxjB 


(33) 


LOADING COEFFICIENT 


By definition the loading coefficient is 


Ap_ 




(34) 


If the body is a cylinder so that its radius does not vary 
with x, the pressure coefficient reduces to 


whore the subscripts L and U refer to the upper and lower 
surfaces of the ah-plane, respectively. It is immediately 
apparent from an inspection of equations (34) and (30) that 
the loading is not affected by A{x ) . Hence, the lift, pitching 
moment, rolling moment, and induced drag for slender 
shapes having a continuous variation of S''{x) can all be 
expressed entirely in terms of ^ 2 (x,*y, 2 ). 

Consider again the type of shapes represented in figure 
3(c) and let there be no discontinuities in S"{x). The 
velocity potential (p 2 for such a class of configurations has 
been expressed in equation (31) as the sum of three poten- 
tials : one due to angle of attack, one due to the thickness of 
the ^ving, and one due to the thickness of the body. It is 
now useful to remark that (pa has odd symmetiy with refer- 
ence to the 2=0 plane and and <po have even symmetry. 


Placing equations (32a) and (32b) into equation (34) and 
using these properties, one finds 


and 


(T Xinr Uo ^ G? dy )Lo 

\ 2 CIo L 


'^<pa J 1_ dR ^ ^ by?g by?c^"[ 


ba; ‘ dx "\bd beJjB 


(35b) 


where A indicates the difference between a quantity on 
vertically opposed pomts of the upper and lower surface of 
the airplane. 

It is apparent from the last two equations that, in general, 
the angle-of-attack and thickness solutions have a coup ling 
effect on the loading coefficient and thei*efore their contribu- 
tion to the load distribution cannot be treated separately. 
It is also important to notice the two special cases in which 
the coupling effects vanish; namely, a body of revolution 
without wings, and an airplane with a cylindrical body 
between the foremost and rearmost extent of the wing. In 
the former case the term b<pcfbd is zero and in the latter 
dRjdx is zero. In both these cases the equation for the 
loading coefficient is 


— =2aA =w A 

2 ox Uo ox 

THE TOTAL LIFT 


(39) 


Total lift can be obtained, of course, by integrating the 
loading coefficient over the aircraft surface. A much simpler 
way of finding the lift, however, can be derived from a 
momentum balance. Thus, by momentum considerations 

it is possible to show that the vectorial force jP on a body 
inside a control surface S is given by the surface integral 

^=- J J p (y-Vo) [v- czs] 

where vector notation is used, the 0 subscript indicates free- 
stream conditions, p and p are the local static pressure and 

density, and V is the local velocity vector. Let the surface 
be a cylinder of infinite radius and two yz planes closing the 
cylinder be located infinitely far ahead of and behind the 
airplane. Then'thelift force is given to the lowest order by 




[pw(?7o+'u)]i=a. dydz 


which, reduces to 


J co 

j {w)x-a. dydz 
-*00 J — oo 

This can be simplified ^ce v)=b<plbz and (A^)x- „ is the same 
as the jump in the potential evaluated at the airplane trailing 
edge. Thus the expression for lift becomes 


=PoJ7of 

J tjy 


(A(p)tb dy 


(37) 
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Equation (37) applies to all slender shapes. In special 
cases represented by figure 3 (c), thickness effects always have 
even symmetry with respect to the 2=0 plane, and it follows 
that the total lift and the vortex distribution in the wake of 
such configurations are affected only by the part of the poten- 
tial due to angle of attack, even though the detailed load 
distribution depends upon both thickness and an^e-of-attack 
solutions. 

EXAMPLES 

PBESSURB ON A THIANQULAR WING WITH ELLIPTIC CROSS SECTION 

It is of interest to calculate, by equation (24), the pressure 
on nonlifting wings of triangular plan form and elliptic cross 
section flying at supersonic speeds, since examples of this 
type have been solved without restriction to slender-wing 
theory. It is proposed, therefore, to study two cases given 
first by Squire (ref. 8) and then to compare the analytical 
results. 

Let the wing be placed at zero angle of attack in a super- 
sonic free stream of Mach number Mo. Consider first the 
thickness distribution for which the ordinate of the upper 
surface is 

where Co is the root chord, i is wing thickness at x=Coj and m 
is the tangent of the semiapex angle of the plan form. The 
flow is supersonic, so it is unnecessary to consider closure. 

It follows from equation (38) that the elliptic section in the 
plane r=xi has major and minor semiaxes equal to mxi 
and txil2co respectively. The cross-sectional area and the 
surface slope are, therefore, 

S{x) =Tdm7^l2ci), bhfbx=tmxl2co V 2 /^ 

Since attention is confined to symmetric nonlifting wings, 
the boimdary conditions are planar and are expressed by 
equation (26) for a=0. Further, the solution is given in 
terms of these boimdary values by equation (23) wherein 
(&9/d7i),«o becomes Uo(^hfdx) and (A^),_o is zero by sym- 
metry. Let us first study the flow on the surface of the wing. ' 
Then, although S'^(x) is discontinuous at the origin, we can 
still write the equation for <p(x, y, 0) in the form (see the 
discussion succeeding equation (24)) 


Since 

1 tmxln\y-y,\dy,_mtx,^mx. . , 

’rJ_m*2coV wV— 2co 2 

the expression for perturbation potential becomes 

ly| < mx 

From equation (32a), wherein dBjdx is, of course, zero 
since there is no body, the pressure coefficient on the wing 
is 


Hence, the pressure distribution on the wing is uniform. 
Analysis not limited by the assumption of slenderness yields 
for pressure coefficient on the wing 

n K—E 

Co Vi— 

where K and E are complete elliptic integrals with modulus 
Since for values of the modulus near one the 
asymptotic relations 


apply, the pressure coefficient in slender-wing theory is seen 
to be a first-order approximation. 

If one is interested in the pressure coefficient in the 2=0 
plane but off (although still in the vicinity of) the iving, the 
discontinuity in S"(x) at the origin should bo considered. 
In such a case, since AjS"( 0) is Td/mjco and is zero, 

equation (24) becomes 




S/I '”[■ 


Xi+V(a;— Xi)^— <3V 


J dxi 


Since 


1 tmxin\y-yMyi ^ ^^ |y|+Vy"-mV . 

one can show the expression for the pressure coefficient off, 
but in the plane of, the wing is 


a 


\y\ 

Co jyj'lf— 771^7? 


-In 


2(a;+V^^n^) . 


\y\>mx 

(40) 


Notice that the pressure has a square-root singularity along 
the wing leading edges, a result consistent with the exact 
lineafized-theoiy solution- The pressure along the Mach 
cone from the wing apex, which is zero according to linearized 
theory, is not zeio according to equation (40) but its magni- 
tude is of the order 05m) ^ which, in keeping with the assump- 
tions of slender-airplone theory, is negligible. 

Squire has also considered the wing with ordinates given by 


A(r,2/)= 



VmV— p 


(41) 


The lateral section is again elliptic, 'with semimajor and 
semiminor axes equal to mx and te®/2co®. Cross-sectional 
area and surface slope are, respectively, 


ci( bh__ i 2mV— p 

2co^ ^ c>x 2mcoVm^'?--“p 

By direct integration it can be shown that on the 'wing 



(39) 


(42) 
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Analysis not limited to slender wings yields the expression 

and again the results are in agreement if higher-order terms 
in /3m ore neglected. 

A study of generalized conical flew flelds in linearized 
supersonic theory reveals that the linear pressme distribu- 
tion in the above problem can also be obtained on a wing 
wth thickness specified by the relation 

Ux,y)=^ cosh-‘ — 

where fc is a constant that can be related to the maximum 
thickness ratio (attained along the line of the 

^ving. The cross-sectional area and surface slope are, 
respectively, 


3 ^ 5a; m^J1n^x^—y^ 

and pressure coofldcient on the wing is 

0,=2kmx(ln^-2^ 

The latter expression agrees to the first order in |3m with the 
general Imoorized solution for such a wing presented in 
reference 9. Slender-wing theory thtis retains the property 
of the more general linear theory in that a given pressure 
distribution does not necessarily yield a unique thickness 
distribution. 

SUPERSONIC DRAG OF WINGS AT ZERO INCIDENCE 

The general expression for the supersonic drag of a slender 
aerodynamic shape has been derived by Word (ref. 3) through 
the use of momentum methods. It is also possible to obtain 
these results by direct integration of the product of pressure 
and surface slope over the specified surface; the analysis, 
however, requires rather careful attention to orders of inte- 
gration, when planar problems are involved. Consider, for 
example, the drag of a wing at zero incidence and with a 
specified thickness distribution s=±A(a;, y). The drag of 
the wing is expressible in the form 

D=D,+2z jj ^dxdy (43) 

whore the first term includes possible contributions to the 
drag that result from a finite leading-edge radixis of curva- 
ture. Prom reference 10, this drag per unit of span is, in 
slender-wing theory. 



(44) 


whore r, is the radius of curvature normal to the wing 
leading edge and s is the local semispan. If the ordinate of 
the wing, in the vicinity of the leading edge, is 


Zu=f(8,y)-Js^ 


equation (44) becomes 


dPe 

dy 


^2 


2 



Assuming the wing is pointed and the only discontinuity 
in the interval 0<a;<Z occurs at the origin, the potential*of 
the wmg, evaluated in the plane of the wing, is given by 






(45) 

and, since pressure coefficient in the planar case is directly 
proportional to the streamwise gradient of <p, the contribu- 
tion of each of the terms on the right-hand side of equation 
(45) can be calculated separately in equation (43). The 
second and third terms offer no difficulty but simplification 
of the expression resulting from the first term necessitates 
an inversion of order of integration and, if the leading edge 
has a finite radius of curvature, such an inversion cannot be 
carried out in the conventional manner. However, a method 
by means of which such an inversion can be carried out is 
presented in reference 11. Thus, set 

and 

where ^ refers to the “finite part” of the integral ’’ and the 


notation J dyj dyi signifies that the yt int^ation must be 
performed first. Then if 


z»(a:, !/)=/(«, 2/) 


it can be shown that 

h 




Detailed analysis reveals that the residual term (i. e., the 
value of h — is) yields a drag component that is equal in 
magnitude but opposite in sign to D,. 

The final expression for the drag of the wing is then 


^ In |-l-i S'(l) f' In Q-Xi) dx,- 

Z'K p TT Jo 

A JJ S"(x)dx JJ S"(Xi) ln\x-xx\dx^ 


+ 


(46) 


As a particular example, consider the wing-like surface of 
triangular plan form (ref. 8) which results from a combina- 


J For a definition of the finite-part IntegratIoh'te<dinIque as used here, see reference 11. 
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tioR of the surfaces specified in equations (88) and (41) and 
has ordinates given by the expression 

Ma;, (co-x)VwV^ (47) 

This wing has rounded leading edges and a finite trailing- 
edge angle, and from equation (46) its drag coeflScient based 
on wing area is found to be 

C7^=-2m,r(lJ(l+Zn^) ' (48) 

■ It IS apparent from equation (46) that, for the type of wing 
considered, wing drag varies with Mach number so long as 
the streamwise gradient of area is finite at the rear of the 
^ving; conversely, there is no dependence on Mach number 
when the gradient of area vanishes there. For example, a 
wing with an elliptic plan form and biconvex sections satisfies 
the latter condition, and its drag coefficient based on wiog 
area is 



where t is total maximum thickness, a is the semiaxis of the 
elliptic plan form in the stream direction, and b is the semi- 
axis measured normal to the stream direction. 

A comparison between the values of given by slender- 
airplane theoiy for the Squire wing (eq. (48)) and the elliptic 
lens (eq. (49j) and the exact thin-airfoil-theory values ® for 
the same wings is shown in figure 4. 


teristics. An essential feature of Jones' results involves 
the use of combined flow flelds that are obtained by super- 
imposing the disturbance fields in forward and reversed 
motions. So long as the governing equations of flow are 
linear, it is possible to establish reciprocity relations between 
the induced fields of arbitrarily situated sources and doublets 
in combined flow fields. Conditions for mmimiim, drag 
under imposed restrictions are then expressed in terms of the 
pressure induced in the superimposed fields. For example, 
it is found that if the thickness distribution for a symmetrical 
nonlifting wing yields a specified volume, then drag is a 
minimum if the thickness is distributed in such a Avay that 
the pressure gradient in the combined field remains constant 
over the plan form of the wing. The application of this 
condition to a slender wing having a continuous value of 
S"{x) eveiywhere except at the origin is simple, since from 
equation (24) the perturbation potential in the two direc- 
tions of flow can be written explicitly. In the plane of the 
wing, the forward flow yields 

¥’/(a;, 2 /, 0 )=^J_ <S"'(a;i) dzi 

and the reverse flow jaelds 

?>r(a:,y,0)=A S"{xi) dxt 

where I is the streamwise length of the wing and b is local 
senuspan. Since 

b(Pf_ b<pr 

bn bn 



Figubb 4. — Comparison of slender-airplane theory with exact thin- 
airfoil theory results for t wo wings in supersonic flight . 


the perburbation potential 4>z in the combined field is 

0)=^£s"(.x,) In dx, 

and it follows directly that if S'(0)=S'(r)=Q, pressure 
coefficient in the combined field is 


0 1 f-SIM 

^ ^'Jo X—Xi 


dxi 


( 60 ) 


In the case of thickness distribution vnth given volume, Op^ 
is a linear function in x and equation (50) is precisely the 
same integral equation that arises in the deterrnination of 
thickness distribution with given volume for a slender body 
of revolution in superonic flight (refs. 14, 16, and 16). The 
same chordwise distribution of area therefore exists for 
wings and bodies of revolution under the given conditions. 
Ames Aeronautical Laboratory 
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